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CHAPTER  1  INTRODUCTION  AND  SUMMARY 


A  portion  of  the  studies  of  the  biological  effects  of  blast  from 
nuclear  explosions  has  been  concerned  with  the  translational  effects 
of  blast  waves  for  objects  as  small  as  a  10-mg  stone  and  as  large 
as  a  168-lb  man.  Computed  results  from  theoretical  studies^-* 
when  compared  to  field  data  for  near-ideal  blast  waves  from  nu¬ 
clear  explosions^,  4,  5  have  demonstrated  that  the  motion  of  experi¬ 
mental  objects  can  be  satisfactorily  predicted  for  free-field  con¬ 
ditions  or  for  window  glass  in  houses. 

This  report  presents  for  high  explosives  (surface  burst)  the 
results  of  a  similar  theoretical  study  —  specifically,  computed 
velocity,  displacement,  and  acceleration  as  functions  of  time  for 
a  variety  of  objects  exposed  to  blast  waves  with  12  maximum  over¬ 
pressures  ranging  from  1  to  20  atm.  Although  all  computations 
were  made  for  1  ton  of  high  explosives,  the  results  may  be  readily 
scaled  to  lower  or  higher  yields.  The  translated  objects,  or  mis¬ 
siles,  are  identified  in  this  study  by  their  acceleration  coefficients* 
which  range  from  0.  01  to  6.  0  ft^/lb.  ** 


*  Acceleration  coefficient  is  defined  for  an  object  as  its  area 
presented  to  the  wind  times  its  drag  coefficient  divided  by  its 
mass.  See  Ref.  2. 

**This  range  in  acceleration  coefficients  is  for  1  ton  of  high 
explosives.  Because  of  scaling  laws  (see  Sect.  2.  1),  different 
ranges  would  apply  to  other  yields. 
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CHAPTER  2  ANALYTICAL  PROCEDURES 
2.  1  THE  MODEL 


The  computational  model  used  in  this  work  was  reported  in 
Ref.  1  and  will  not  be  described  in  detail.  In  the  previous  study 
tabular  values  of  computed  velocity,  displacement,  and  accel¬ 
eration  as  functions  of  time  were  presented  as  nondimensional 
quantities  for  missiles  produced  by  nuclear  blast  waves.  How¬ 
ever,  to  make  interpretation  of  the  results  easier,  the  compu¬ 
tations  in  the  present  study  for  missiles  produced  by  high- 
explosive  blast  waves  were  made  in  dimensional  form  for  a  yield 
of  1  ton,  and  the  results  are  presented  graphically.  Dimensional 
analysis  derived  previously^  make  it  possible  to  apply  the  results 
to  explosions  of  lower  or  higher  yields  within  the  limits  of  weapons' 
scaling. 

In  deriving  the  computational  model  (see  Ref.  1),  the  following 
assumptions  were  made:  (1)  Friction  between  the  missile  (trans¬ 
lated  object)  and  the  surface  was  negligible.  (2)  The  effect  of 
gravity  on  the  horizontal  velocity  of  the  missile  could  be  neglected. 
(3)  Only  the  winds  associated  with  the  blast  wave  contributed  to 
missile  translation.  (4)  The  acceleration  coefficient  of  the  mis¬ 
sile  could  be  assumed  to  be  constant.  (5)  The  blast  wave  does  not 
decay  appreciably  while  passing  the  missile. 

Allowance  was  made  for  the  object's  being  exposed  to  the  blast 
winds  for  a  time  dependent  on  the  relative  velocity  of  the  missile 
and  the  blast  wave.  This  effect  was  particularly  important  for  the 
missiles  with  the  higher  acceleration  coefficients  which  are  pro¬ 
pelled  to  relatively  high  velocities. 

Numerical  solutions  of  the  model  were  determined  by  stepwise 
integration  of  the  model  equations.  1  For  missile  velocity,  the 
following  was  used 

±  Av  =  e  +  f  -J e^  +  f^+  2  f  g  (1) 

where 


Av  is  positive  if  u  >  v  , 
r  o  o 

Av  is  negative  if  uq  <  vq, 

Av  =  change  in  missile  velocity  during  time  step  At, 


e  =  x  -  v  , 
—  o 
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f  =  a  q  (u  -  v  )  x  At/u  , 
g  =  x  -  u, 

uq  =  wind  velocity  at  beginning  of  At, 
vq  =  missile  velocity  at  beginning  of  At, 
u  =  average  wind  velocity  during  At, 

x  =  average  velocity  of  propagation  of  blast  wave  during  At, 

—  .  f~2  2 

approximated  by  x  =  0.  6  u  +  /  cq  +  .  36  _u 

cq  =  velocity  of  sound  in  the  undisturbed  air 

a  =  acceleration  coefficient  =  s  C,/m, 

d 

s  =  area  presented  to  wind  by  missile, 

=  drag  coefficient  of  missile, 

m  =  mass  of  missile,  and 

q  =  average  dynamic  pressure  during  At. 

Incremental  distance,  Ad,  was  computed  by  the  following: 

Ad  =  (v  +  Av/2)  At  x/  (x  -  vq  -  Av/ 2)  (2) 

Missile  acceleration,  a,  was  determined  from  the  following 
(integration  being  unnecessary): 


a  =  ‘ar  ~  9. a  (u  ■  v)2  /  u2  (3) 

where  q,  u,  and  v  are  dynamic  pressure,  wind  velocity,  and  missile 
velocity  at  any  time  t.  (Note  that  an  acceleration  numeric  is  used  in 
Chap.  3:  A  =  a/g  where  g  is  the  acceleration  of  gravity.  ) 

Equations  (1)  and  (2)  were  integrated  in  a  stepwise  fashion  from 
the  arrival  time  of  the  blast  wave  (t  =  0)  to  the  time  of  zero  wind 
velocity  (t  =tj).  Because  of  rapid  changes  in  missile  and  wind 
velocities  shortly  after  the  arrival  of  the  blast  wave,  it  was  neces¬ 
sary  to  use  smaller  time  steps  during  the  early  times  than  during 
later  times.  The  following  arbitrary  scheme  was  used  to  determine 
the  variable  time  step.  The  first  step  was  always  0.  001  tj.  The 
remaining  time  (tj  -  0.  001  tj)  was  divided  into  99  log  intervals. 

The  first  85  of  the  log  intervals  were  used  as  ever-increasing  time 
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steps.  The  85th  log  interval  was  then  used  as  a  constant  time  step 
until  time  tj  was  reached. 

For  convenience,  the  scaling  laws  for  translational  studies 
derived  in  Chaps.  2  and  5  of  Ref.  1  will  be  restated  using  the  termi¬ 
nology  of  the  present  report.  The  subscript  "1"  is  used  to  denote 
parameters  applicable  to  a  yield  of  1  ton  of  high  explosives,  to  an 
ambient  speed  of  sound  of  1117  ft/ sec,  and  to  an  ambient  pressure 
of  14.7  psi.  The  parameters  without  subscripts  are  applicable  to  a 
yield  of  W  tons,  to  an  ambient  speed  of  sound  of  cD  ft/ sec,  and  to 
an  ambient  pressure  of  pQ  psi.  Thus,  the  results  in  Chap.  3  can 
be  scaled  as  follows: 

v  =  v,  (c  /  1 1 1 7 ) 

1  o 

d  =  dt  (14.  7  W/po)1/3 
A  =  (c  /  1 1 1 7) 2  (1/W)1/3  (p  /  14.  7)i/3 
t  =  t1  (1 1  17/cq)  (14.  7  W/po)1/3 
a  =  (14.7/pq)2/3  (co/1117)2  (1/W)1/3 

where 


v  =  missile  velocity  in  ft/ sec, 

cq  =  ambient  speed  of  sound  in  ft/ sec, 

d  =  distance  of  missile  travel  in  ft, 

W  =  yield  of  high  explosives  in  tons, 
pQ  =  ambient  pressure  in  psi, 

A  =  acceleration  of  missile  (A  =  a/32.  2)  in  gravity  units, 
t  =  time  after  arrival  of  blast  wave  in  msec, 

2 

a  =  acceleration  coefficient  (c  =  s  C^/m)  in  ft  /lb. 

2.2  BLAST- WAVE  PARAMETERS 
2.  2.  1  General 

The  solution  of  the  translation  model  described  in  the  last 
section  requires  that  dynamic  pressure  (q)  and  wind  velocity  (u)  be 
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defined  as  a  function  of  time.  Clear  and  explicit  presentations  of 
these  quantities  for  high  explosives  were  not  found  in  the  literature. 
Consequently,  most  of  the  blast  material  used  in  this  report  was 
taken  from  a  numerical  study  by  Brode^  and  from  experimental 
results  reported  by  Goodman.  3  This  material  along  with  Rankine- 
Hugoniot  equations^  allowed  computation  of  the  needed  blast  para¬ 
meters.  The  overpressure-time  relation,  needed  to  determine 
wind  velocity  as  a  function  of  time,  will  be  treated  first. 

2.  2.  2  Overpressure  vs.  Time 


The  overpressure  information  was  obtained  from  experi¬ 
mental  results.  3  Duration  and  overpressure  impulse  scaled  to 
1  ton  of  high  explosives  are  recorded  in  Table  2.  1  for  the  overpres¬ 
sure  values  of  interest  in  this  study.  Overpressure  as  a  function 
of  time  was  not  defined  in  Ref.  3,  but  in  the  present  report  it  was 
assumed  to  be  of  the  form 

P  =  P  (1  -  T)e-nT  (4) 

s 

where 

P  =  overpressure  in  atm, 


Pg  =  maximum  or  shock  overpressure  in  atm, 

T  =  t/t+, 

P 

t  =  time  after  arrival  of  the  blast  wave  in  msec, 


tp  =  duration  of  positive  overpressure  in  msec,  and 
n  =  a  constant  for  a  given  value  of  P  . 

Since  n  in  the  above  equation  determines  the  shape  of  the 

rt+ 

p 

P-t  curve,  it  also  determines  the  impulse,  QJ  P  dt,  for  particular 

values  of  Ps  and  t£.  Thus,  integrating  Eq.  (4)  gives  the  following 
impulse,  Ijij,  relation: 

I+  =  P  t+  (e  n  +  n  -  i)/n^  (5) 

P  s  p 

The  values  of  n  listed  in  Table  2.  1  were  found  using  Eq.  (5). 
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Table  2.  1 


Parameters  for  Determining  Overpres sure -time  Functions* 
(1  Ton  of  High  Explosives,  Surface  Burst) 


P  ,  atm 
s 

,  + 

t  ,  msec 

P 

i"*",  atm  msec 

P 

n 

1.  0 

16.  5 

6.  48 

.  960 

1.  5 

14.  4 

7.  55 

1.  32 

2.  0 

12.  9 

8.  30 

1. 47 

2.  5 

11.  7 

8.  91 

1. 65 

3.  0 

10.  9 

9.  46 

1.  81 

4.  0 

9.  70 

10.  3 

2.  15 

5.  0 

8.  87 

10.  9 

2.  42 

6.  0 

8.  40 

11.4 

2.  80 

8.  0 

7.  99 

12.  2 

3.  63 

10.  0 

7.  73 

12.  7 

4.  40 

15.  0 

6.  27 

13.  6 

5.  19 

20.  0 

3.  00 

14.  2 

2.  65 

*P  =  P  (1 
s 

rpv  -nT 
-  T)  e 

where  T  =  t/t+ 

P 
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2.  2.  3  Dynamic  Pressure  vs.  Time 


Dynamic  pressure  in  atmospheres  at  the  shock  front,  Qs> 
was  computed  using  the  Rankine-Hugoniot  relation  reported  in  Ref.  4: 

2 

r*  “ 


Qs  = 


2.  5  P 

c 

7  +  P 


Values  of  Qs  corresponding  to  Ps  values  used  in  this  study  are 
listed  in  the  second  column  of  Table  2.  2.  Durations,  tj,  and  im¬ 
pulses,  ij,  for  dynamic  pressure  in  the  same  table  were  obtained 
by  scaling  the  results  of  a  numerical  study^  to  a  yield  of  1  ton 
(surface  burst).  A  scaling  factor  was  applied  to  the  numerical  re¬ 
sults  in  order  to  make  the  overpressure  durations  consistent  with 
those  found  experimentally.  ^ 


A  procedure  similar  to  that  described  in  Sect.  2.  2.  2  was 
used  to  determine  dynamic  pressure  vs.  time.  Values  recorded 
in  Table  2.  2  for  r  were  obtained  from 

I+  =  Q  t+  (e"r  +  r  -  l)/r2  (7) 

u  s  u 

where 


Qg  =  maximum  or  shock  dynamic  pressure  in  atm, 

t^  =  duration  of  positive  dynamic  pressure  in  msec,  and 

r  =  a  constant  for  a  particular  value  of  P  or  Q  . 

s  s 

Dynamic  pressure  as  a  function  of  time  could  then  be  found 

using 

Q  =  Q  (1  -  T)  e'rT 
s 

where  T  =  t/t*. 

u 

2.  2.  4  Wind  Velocity  vs.  Time 


By  definition, 

q  =  1/2  p  u2  (8) 

where  q  =  dynamic  pressure, 
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Table  2.  2 


Parameters  for  Determining  Dynamic  Pressure 
as  a  Function  of  Time* 

(1  Ton  of  High  Explosives,  Surface  Burst) 


P  ,  atm 
s 

Q  ,  atm** 
s 

,+ 

t  ,  msec 
u 

I+,  atm  msec 
u 

r 

1.  0 

.  3125 

23.  4 

1.  56 

3.  33 

1.  5 

.  6618 

20.  4 

2.  14 

5.  10 

2.  0 

1.111 

18.  4 

2.  72 

6.  32 

2.  5 

1. 645 

17.  1 

3.  26 

7.  47 

3.  0 

2.  250 

16.  1 

3.  80 

8.  39 

4.  0 

3.  6 36 

14..4 

4.  80 

9.  80 

5.  0 

5.  208 

13.  1 

5.  77 

10.  7 

6.  0 

6.  923 

11.9 

6.  65 

11.  3 

8.  0 

10. 667 

9.  85 

8.  68 

11.0 

10.  0 

14. 706 

8.  17 

10.  81 

10.  0 

15.  0 

25.  568 

5.  06 

16.  43 

6.  69 

20.  0 

35. 556 

3.  00 

22.  19 

3.  46 

*Q  = 

Qs  (1  -  T)  e"rT 

where 

T  =  t/t+ 
u 

? 

2.  5  P 

g 

** Computed  from  Q  =  - — — — — 

s  7  +  F„ 
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p  =  air  density,  and 
u  =  wind  velocity. 


Wind  velocity  could  thus  be  determined  after  the  dynamic  pressure 
and  the  air  density  were  evaluated.  The  air  density  across  the  . 
shock,  p  ,  was  found  using  one  of  the  Rankine-Hugoniot  equations: 

Ps  =  PQ  (?  +  6  Ps)  / ( 7  +  Ps)  (9) 

where  pQ  is  the  ambient  air  density. 


Changes  in  air  density  after  the  passage  of  the  shock  were  assumed 
to  be  adiabatic.  Thus, 


,P  +  1  x  1/1.  4 
+  V 


s 

where  p  =  air  density  when  the  overpressure  is  P  atm. 


(10) 
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CHAPTER  3  RESULTS 
3.  1  GENERAL 


Computed  results  were  obtained  using  the  translation  model 
described  in  Sect.  2.  1  and  a  digital  computer  with  an  incre¬ 
mental  plotter  to  graph  the  output  data.  The  system  for  deter¬ 
mining  time  steps  (see  Sect.  2.  1)  resulted  in  106  to  112  steps 
for  each  numerical  integration  from  the  arrival  of  the  blast  wave 
(t  =  0)  to  the  time  of  zero  wind  (t  =  t^). 

As  previously  noted,  all  solutions  were  made  for  1  ton  of 
high  explosives  burst  at  the  surface.  The  12  different  blast  waves 
used  are  identified  in  terms  of  overpressure:  Ps  =  ps/p  0  ,  the 
ratio  of  shock  overpressure  to  local  ambient  pressure  (not 
necessarily  the  sea-level  value  of  14.  7  psi),  sometimes  called 
excess  pressure  ratio.  Values  of  Ps  used  were  1.  0,  1.  5,  2.  0, 
2.5,  3.0,  4.0,  5.0,  6.0,  8.0,  10.0,  15.0,  and  20.  0  atm. 

For  each  Ps,  numerical  integrations  were  made  for  the 
following  acceleration  coefficients,  ay,  in  the  order  listed:  6,  3, 

2,  1,  0.6,  0.3,  0.2,  0.1,  0.06,  0.03,  0.02,  0.0ift2/lb.  If  the 
maximum  velocity  computed  for  any  acceleration  coefficient  was 
less  than  10  ft/ sec,  the  computations  were  halted  for  that  over¬ 
pressure. 

According  to  the  translational  model  used,  the  behavior  of 
an  object  is  determined  by  its  acceleration  coefficient,  all  other 
factors  being  constant.  To  aid  interpretation  of  the  computed  re¬ 
sults,  a  list  of  acceleration  coefficients  obtained  from  Ref.  1 
for  a  variety  of  objects  is  reproduced  in  Table  3.  1.  A  more  com¬ 
plete  source  of  acceleration-coefficient  information  can  be  found 
in  Ref.  2. 

3.  2  VELOCITY  VS.  TIME  AND  DISTANCE  VS.  TIME 


Computed  velocity  and  distance  as  functions  of  time  are  pre¬ 
sented  for  a  maximum  overpressure*  of  1  atm  in  Figs.  3.  1  and 
3.  2,  respectively,  for  seven  acceleration  coefficients.  Machine 
plots  for  these  and  succeeding  figures  connected  with  straight 
lines  every  other  computed  point  for  the  first  86  time  steps.  All 
of  the  remaining  time  steps  were  plotted.  Each  of  the  curves  end 
at  the  midpoint  of  the  last  time  step  before  the  dynamic  pressure 


^Defined  in  Sect.  3.  1. 
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Table  3.  1 


Typical  Acceleration  Coefficients,*  a  =  sC^/m 
where  s  is  the  area  presented  to  wind  by  missile, 
is  the  drag  coefficient  of  the  missile, 
and  m  is  the  mass  of  the  missile 


1 68-lb  man: 

c,  ft2 /lb 

Standing  facing  wind 

0.  052 

Standing  sidewise  to  wind 

0.  022 

Couching  facing  wind 

0.  021 

Couching  sidewise  to  wind 

0.  017 

Prone  aligned  with  wind 

0.  0063 

Prone  perpendicular  to  wind 

0.  022 

Average  value  for  tumbling  man 
in  straight,  rigid  position 

0.  03 

21-g  mice,  maximum  area  presented  to  wind 

0.  38 

180-g  rats,  maximum  area  presented  to  wind 

0.  19 

530-g  guinea  pigs,  maximum  area  presented  to  wind 

0.  15 

2100-g  rabbits,  maximum  area  presented  to  wind 

Typical  stones: 

0.  079 

0.  1  g 

0.  67 

1.  0  g 

0.  32 

10.0  g 

0.  15 

Window -glass  fragments,  1/8  in.  thick** 

0.  1  g,  all  orientations 

0.  78 

1.  0  g,  edgewise  and  broadside  to  wind 

0.  48-0.  57 

10.  0  g,  edgewise  and  broadside  to  wind 

0.  34-0. 72 

*F rom  Ref.  1 . 

**S ingle -strength  window  glass.  See  Ref.  2  for  data  on  plate 
glas  s . 
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Fig.  3.  2  Distance  vs.  Time  for  P  =1.0  atm 


and  winds  become  negative.  For  high  acceleration  coefficients, 
the  curves  terminate  at  later  times  than  do  those  for  low  ones: 
the  missiles  with  high  coefficients  —  and,  thus,  with  high  veloc¬ 
ities  —  travel  along  with  the  blast  wave  for  longer  times  than  do 
those  with  the  low  ones. 

The  appropriate  equations  for  scaling  the  computed  results 
to  other  yields  and  ambient  pressures  and  speeds  of  sound  are 
presented  on  each  chart.  To  illustrate  scaling  from  1  ton  to 
1000  tons  (1  kt)  and  to  compare  the  results  for  1  kt  with  those 
for  nuclear  blast  waves,  *  consider  the  maximum  velocity  and 
distance  of  travel  at  maximum  velocity  predicted  for  a  1-gm 
stone  when-  p  =  1 .  0  atm,  p  0  =  14.7psi,  cQ=  1117  ft/  sec.  A 
stone  of  1  gm  has  an  acceleration  coefficient  of  0.  32  ft^/lb. 

(See  Table  3.  1.  )  For  a  yield  of  1  ton  the  maximum  predicted 
velocity  for  =  0.  32  ft^/lb  obtained  from  Fig.  3.  1  is  about 
30  ft/ sec  occurring  19  msec  after  the  arrival  of  the  blast  wave. 

By  referring  to  Fig.  3.  2,  the  distance  of  travel  of  19  msec  is 
found  to  be  about  0.  48  ft. 

To  apply  the  computed  data  to  a  yield  of  1000  tons,  it  is  first 
necessary  to  determine  an  equivalent  acceleration  coefficient, 
a\_,  for  a  yield  of  1  ton.  By  using  the  scaling  equation  for  accel¬ 
eration  coefficient  in  Figs.  3.  1  and  3.  2,  a  j  =  0.  32  ft^/lb  x  (1000)*'  ^ 

=  3.  2  ft^/lb.  The  maximum  velocity  predicted  for  this  value  of 
is  about  185  ft/ sec  occurring  15  msec  after  the  arrival  of  the  blast 
wave  (Fig.  3.  1).  The  distance  traveled  for  W  =  1  ton  and  for 

=  3.  2  ft^/lb  at  15  msec  is  2.  4  ft  (Fig.  3.  2).  For  W  =  1000  tons, 
the  distance  is  2.4  ftx  ( 1000)  *  I  ^  =  24  ft  occurring  15  msec  x 
(1000)*'^  =  150  msec  after  the  arrival  of  the  blast  wave.  For  com¬ 
parison,  the  maximum  velocity  and  distance  of  travel  at  maximum 
velocity  computed  for  a  nuclear  blast  wave*  for  the  conditions  stated 
above  for  a  yield  of  1  kt  are  200  ft/ sec  (high  explosives:  185  ft/ sec) 
and  28.  7  ft  (high  explosives:  24  ft).  Similar  comparisons  for  lower 
acceleration  coefficients,  however,  showed  a  larger  discrepancy 
between  the  velocities  and  distances  predicted  for  the  high -explosive 
and  nuclear  blast  waves.  The  explanation  of  this  is  that  the  blast 
criteria  used  for  P  =  1  atm  indicates  a  shorter  duration  of  dynamic 
pressure  for  the  high-explosive  than  for  the  nuclear  blast  wave  for 
the  same  explosive  yield,  maximum  overpressure,  and  ambient 
pressure  and  speed  of  sound.  Duration  effects  on  maximum  velocity 
are  less  pronounced  for  the  high  acceleration  coefficients  than  for  the 
low  ones  since  the  former  reach  maximum  velocity  before  the  blast 
wave  has  decayed  appreciably. 

The  charts  shown  in  Figs.  3.  3  to  3.  24,  similar  to  those  described, 
were  computed  for  maximum  overpressures  of  1.  5,  2.  0,  2.  5,  3.  0, 
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Fig.  3.4  Distance  vs.  Time  for  P  =1.5  atm 


For  other  conditions  use:  =  a  (1117/02  <po/14.  7)2/3  W1/3;  t  =  (1117/cO  (14.7  W/p  J1/3; 

and  v  =  v,  (c  /1117) 


Fig.  3.  5  Velocity  vs.  Time  for  P  =2.0  atm 


Distance  vs.  Time  for  P  -  2.0  atm 


For  other  conditions  use:  a,=  a  (1117/c  )2  (p  /14.7)2/3  wl/3.  t  =  t.  (Ill  7/c  )  (14. 7  W/p  )l/3; 


For  other  conditions  use:  a,=  o  (1117/c  )2  (p  /14.7)2/3  W1/3;  t  =  t,  (1117/c  H14. 7  VV/p  )V3 
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Fig.  3.  8  Distance  vs.  Time  for  P  =2.5  atm 


400 


23 


Velocity  vs.  Time  for  P  =  3.  0  atm 


For  other  conditions  use:  a,  =  a  (1117/c  >2  (p  /14. 7)2/3  W1/3;  t  =  t,  C1117/cH14.7  W/p  )1/3 


0  atm 


For  other  conditions  use:  a,  =  a  (111  7/c  (p  /14.  7)^/3  \yl/3.  t  =  t,  (1117/c  ) (14. 7  W/p 


Fig.  3.  12  Distance  vs.  Time  for  P  =  4.  0  atm 


0  atm 


For  other  conditions  use:  a,  =  a  (1117/c  (p  /14.  7)^/3  \vV3.  t  =  t,  (1117/c  )  (14. 7  VV/p  )l/3. 


0  atm 


For  other  conditions  use:  a,  =  a  (1117/c  F  (p  /14.  7)2/3  yyl/3.  t  =  t,  (11 17/c  )  (14. 7  \V/p  )l/3 
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Fig.  3.  16  Distance  vs.  Time  for  P  =  6.  0  atm 


For  other  conditions  use:  a,=  a  (1117/c  )2  (p  /1 4. 7)2/3  v\'F3;  t  =  t,  (11 17/c  )  (14. 7  W/p  )V3 


Fig.  3.20  Distance  vs.  Time  for  P  =  10.0  atm 


Time  for  P  =  15.  0  atm 


36 


Fig.  3.22  Distance  vs.  Time  for  P  =  15.0  atm 


•  conditions  use:  cu  =  o  (1117/c  (p  /14. 7)^/3  W-/-*;  t  -  t,  (1117/c  ) (14. 7  \V/p 


ie  for  P  =  20,  0  atm 


For  other  conditions  use:  a.  =  c  (11 17/c  >2  (p /14. 7)2/3  W1/3;  t  =  t,  (1117/c  )  (14. 7  W/p  )V3 


4.  0,  5.  0,  6.  0,  8.  0,  10.  0,  15.  0,  and  20  atms.  The  velocity- 
curves  for  ay  -  6.  0,  in  some  cases,  cross  those  for  the  lower 
acceleration  coefficients  (see  Fig.  3.  9,  for  example).  Com¬ 
parison  of  the  velocity  charts  with  the  -  last-wave  parameters 
shows  that  this  phenomenon  occurs  only  for  blast  waves  whose 
dynamic  pressure  decays  relatively  fast  with  time;  i.  e.  ,  blast 
waves  identified  with  the  higher  values  of  "r11  listed  in  Table  2.  2. 

3.  3  ACCELERATION  VS.  TIME 

When  a  blast  wave  first  encounters  a  missile  having  zero 
velocity,  the  maximum  acceleration  experienced  by  the  missile 
is  the  product  of  its  acceleration  coefficient  and  the  maximum 
dynamic  pressure.*  Computed  values  of  maximum  acceleration 
in  g-units  are  presented  in  Table  3.  2  for  the  maximum  over¬ 
pressures  and  for  the  acceleration  coefficients  used  in  this  study. 
After  the  missile  attains  a  finite  velocity,  however,  missile  and 
wind  velocities  also  control  missile  acceleration.*  Thus,  scaling 
procedures  are  not  nece  ssary  to  obtain  maximum  acceleration, 
but  scaling  (as  indicated  on  Figs.  3.  25  to  3.  36)  is  necessary  for 
evaluation  of  accelerations  occurring  after  the  maximum. 

Plots  of  acceleration  vs.  time  are  in  Figs.  3.  25  to  3.  36  for 
the  same  combinations  of  overpressure  and  acceleration  coef¬ 
ficients  for  which  velocity  and  displacement  data  were  presented 
in  the  last  section.  In  order  to  separate  the  curves  appearing  on 
each  chart,  the  plots  were  not  always  made  to  zero  acceleration. 


*This  relation  is  expressed  in  Eq.  (3)  in  Chap.  2: 
a  =  q  a  (u  -  v)^/u^. 
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Table  3.  2 


Maximum  Acceleration  (A,  g-units)  for  12  Acceleration  Coefficients 
(a,  ft“/lb)  and  for  12  Maximum  Overpressures  (P  ,  atm) 


p 

s 

a  -  6  a  -  3 

a  =  2 

a  =  1 

a  =  .  6 

a  =  .  3 

a  =  .  2 

1.  0 

3,970  1,980 

1,  320 

662 

397 

198 

132 

1.  5 

8,410  4,200 

2,  800 

1, 400 

841 

420 

280 

2.  0 

14,  100  7,  060 

4,  700 

2,  350 

1,410 

.  706 

470 

2.  5 

20,900  10,400 

6,  960 

3,  480 

2,  090 

1,  040 

696 

3.  0 

28,  600  14,  300 

9,  530 

4,  760 

2,  860 

1,  430 

953 

4.  0 

46,  200  23,  100 

15, 400 

7,  700 

4,  620 

2,  310 

1,  540 

5.  0 

66,  100  33,  100 

22, 000 

1  1,  000 

6,  610 

3,  310 

2,  200 

6.  0 

87,900  44,000 

29,  300 

14, 700 

8,  790 

4,  400 

2,  930 

8.  0 

135,000  67,700 

45, 200 

22, 600 

13,  500 

6,  770 

4,  520 

10.  0 

187,000  93,400 

62, 300 

31,  100 

18, 700 

9,  340 

6,  230 

15.  0 

325,000  162,000 

108,  000 

54, 100 

32, 500 

16, 200 

10, 800 

20.  0 

452,000  226,000 

151, 000 

75, 300 

45, 200 

22, 600 

15,  100 

P 

S 

a  -  .  1 

a  -  .  06 

a  -  . 

03  a 

.  =  .  02 

a  =  .  01 

1.  0 

66.2 

39.  7 

19. 

8 

13.  2 

6.  62 

1.  5 

140 

84.  1 

42. 

0 

28.  0 

14.  0 

2.  0 

235 

141 

70. 

6 

47.  0 

23.  5 

2.  5 

348 

209 

104 

69.  6 

34.  8 

3.  0 

476 

286 

143 

95.  3 

47.  6 

4.  0 

770 

462 

231 

154 

77.  0 

5.  0 

1,  100 

661 

331 

220 

110 

6.  0 

1, 470 

879 

440 

293 

147 

8.  0 

2,  260 

1,  350 

677 

452 

276 

10.  0 

3,  1  10 

1,  870 

9  34 

623 

311 

15.  0 

5,  410 

3,  250 

1, 620 

1,  080 

541 

20.  0 

7,  530 

4,  520 

2,  260 

1, 

510 

753 
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Fig.  3.  25  Acceleration  vs .  Time 
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Computed  for:  P  =  2.0  atffl;  VV-lton;  p  — 14.7psi;  and  c  - 1117  ft/sec 


ation  vs.  Time  for  P  -2.0  atm 


45 


ation  vs.  Time  for  P  -  3.  0  atm 
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Fig.  3.  32  Acceleration  vs.  Time  for  P  =6.0  atm 
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34  Acceleration  vs.  Time  for  P  =  10.0  atm 
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Fig.  3.35  Acceleration  vs.  Time  for  P  =  15.0  atm 
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Time,  t^,  msec 

Fig.  3.  36  Acceleration  vs.  Time  for  P  =  20.  0  atm 
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CHAPTER  4  DISCUSSION 


4.  i  MODEL  RELIABILITY 


This  report  presents  numerical  predictions  of  the  behavior 
of  objects  set  in  motion  by  high-explosive  blast  waves.  Un¬ 
fortunately,  these  predictions  cannot  be  compared  with  experi¬ 
mental  data.  The  translation  model,  *  however,  has  been 
successfully  used  to  predict  the  results  of  secondary-missile 
experiments  made  at  the  Nevada  Test  Site  with  nuclear -produced 
blast  waves  which  were  near  ideal  (or  classical)  in  character.* 

In  another  experiment  with  anthropomorphic  dummies^’  maximal 
velocity  could  be  successfully  computed  using  an  average  acceler¬ 
ation  coefficient  for  a  tumbling  dummy.  However,  to  duplicate 
more  precisely  the  velocity-distance  measurements,  it  was  necessary 
to  use  an  acceleration  coefficient  which  was  a  function  of  the  orien¬ 
tation  of  the  dummy  during  translation. 

4.  2  COMPARISON  WITH  NUCLEAR  TRANSLATIONAL  EFFECTS 


In  comparison  with  results  computed  for  nuclear  blast  waves, 
those  for  high-explosive  waves  indicate  that  the  overpressure  must 
be  considerably  higher  for  an  object  to  attain  the  same  maximum 
velocity.  This  velocity  occurs,  however,  after  a  much  shorter  dis¬ 
tance  of  translation.  Because  of  the  short  distances  involved,  it 
seems  reasonable,  in  many  cases,  to  assume  that  a  translated  man 
would  not  change  orientation  during  the  accelerative  phase  of  dis¬ 
placement  induced  by  high  explosives;  thus,  a  nonvarying  acceleration 
coefficient  corresponding  to  that  of  his  original  posture  could  be  used. 
For  example,  the  charts  in  Fig.  3.  13  and  3.  14  for  Pg  =  5  atm  and 
W  =  1  ton  show  that  a  standing  person  with  an  acceleration  coefficient 
of  0.  06  ft^/lb  would  attain  a  velocity  of  23  ft/  sec  in  only  0.  1  ft  of 
travel. 


A  comparison  was  made  in  Sect.  3.  2  between  the  velocities 
predicted  using  nuclear  and  high-explosive  blast  data  evaluated  for 
the  conditions:  P_  =  1  atm,  W  =  1  kt,  I?_  =  14.  7  psi,  c  =1117  ft/sec. 


X 


t 


*In  these  experiments,  reported  in  Ref.  2,  the  velocities 
were  measured  for  stones  and  spheres  in  open  areas  and  for  glass 
fragments  from  windows  facing  the  oncoming  blast  wave.  The  blast 
wave  entering  the  houses  through  the  windows  was  modified;  how¬ 
ever,  if  it  was  assumed  to  have  a  maximum  overpressure  equal  to 
the  reflected  value  of  normal  incidence,  the  maximum  fragment 
velocities  could  be  predicted. 


-54- 


¥ 


For  high  acceleration  coefficients  the  velocities  were  in  reason¬ 
able  agreement,  but  not  for  low  coefficients.  This  discrepancy 
results  from  differences  in  the  blast  parameters  specified  for 
the  two  types  of  blast  waves  and  serves  to  emphasize  the  im¬ 
portance  of  input  data  in  determining  the  computed  results.  One 
method  for  verifying  the  input  data,  as  well  as  the  model  used  in 
the  present  study,  is  to  perform  field  experiments  similar  to 
those  done  with  nuclear  explosions.  ^ 

4.  3  BIOMEDICAL  INTERESTS 


4.  3.  1  General 


Those  interested  in  the  relation  between  environmental 
medicine  and  weapons  effects  recognize  that  any  reasonably 
complete  understanding  of  the  many  problems  envolved  requires 
information  in  the  physical,  biophysical  and  biomedical  areas. 

In  this  regard,  a  conceptual  guide  for  analytical  procedures  and 
research  planning  is  essential;  indeed  such  has  been  proposed^ 
wherein  five  problem  areas  were  defined  to  elucidate  the  kinds 
of  data  needed  to  establish  a  quantitative  fabric  that  would  allow 
the  source  of  an  environmental  variation  to  be  "tied"  to  hazards 
assessment. 

The  five  problem  areas,  plus  another  concerned  with  bio  ¬ 
medical  tasks,  are  listed  in  Table  4.  1.  The  first  three  —  encom¬ 
passing  "free-field"  scaling,  "geometric"  scaling  and  secondary 
events  —  represent  ground  that  must  be  11  spaded1'  mostly  by  those 
qualified  in  the  physical  sciences  if  understanding  of  the  environ  ■ 
mental  variations  that  can  occur  at  potentially  populated  locations 
is  to  be  forthcoming. 

Contemplation  of  the  remaining  three  problem  areas  make  it 
apparent  that  hazards  assessment  requires  knowledge  of  biologic 
response  and  the  etiologic  mechanisms  involved.  Such  knowledge, 
in  turn,  touches  biomedical  tasks  such  as  therapy,  rehabilitation 
and  all  possible  means  for  minimizing  casualties  through  what¬ 
ever  protective  measures  might  prove  effective  and  feasible.  It 
is  here  that  personnel  qualified  in  biophysics,  biology  and  medicine 
can  contribute. 

4.  3.  2  The  Translational  Problem 

Missiles 

Since  experience  has  shown  that  blast-induced  environmental 
variations  which  are  potentially  hazardous  include  the  translation 
of  both  animate  and  inanimate  objects,  applicable  and  definitive 
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Table  4.  1 


Problem  Areas  Relevant  to  Biologic 
Effects  of  Nuclear  Weapons 


Source 

Design 

Yield 

Burst  conditions 

Range 

W  eather 

"Free  -field" 

scaling* 

Attenuation  and 

Augmentation 

Modification  of  "free- 

field"  phenomena  by 

geometric  conditions 

of  exposure 

"Geometric" 

scaling 

Physical 

Interaction 

Energy  transfer  to: 

Physical  objects  and 
biological  material 

Secondary 

events 

Biophysical 

Interaction 

Energy  dissipation  by 
or  within  biologic 

targets 

Etiologic 

mechanisms 

Biologic 

Response 

Major  medical  syndromes 

Isolated  individual  effects 

and  combined  injury 

Hazard 

assessment 

Biomedical 

Tasks 

Therapeutic  and  prophy¬ 
lactic  measures 

Casualty  care 

Rehabilitation 

Protective  pro¬ 
cedures 

*See  Fig.  4.  1  which  shows  the  maximal  values  of 
overpressure  as  a  function  of  range  from  a  1  ton  surface  burst  of 
high  explosives  at  sea  level.  The  chart  is  useful  since  it  allows 
one  to  determine  ranges  for  the  overpressures  mentioned  in 


Sect.  3.1. 
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data  are  required  for  any  comprehensive  analysis.  For  example, 
among  the  factors  that  contribute  to  the  casualty  potential  of  blast- 
energized  missiles  are  the  velocity  and  angle  of  impact;  the  mass 
density,  shape  and  character  of  the  debris;  and  the  area  of  the 
body  receiving  penetrating  and/or  nonpenetrating  wounds. 

Displacement 


Similarly,  the  potential  for  injury  as  a  consequence  of  gross 
displacement  of  biological  targets  may  be  due  to  accelerative 
and/or  decelerative  loading.  The  former  depends  at  least  upon 
the  magnitude  of  the  forces  vs.  time  which  initiate  displacement 
and  upon  the  initial  and  subsequent  orientations  of  the  biologic 
target.  The  latter  depends  mostly  upon  the  velocity  at  which 
deceleration  occurs,  the  character  of  the  decelerating  surface 
and  the  area  or  areas  of  the  body  involved  whether  impacting 
with  a  solid  object  or  tumbling  over  some  near  horizontal  sur¬ 
face  transpires. 

4.  3.  3  Present  Study 

The  previous  paragraphs  help  to  place  in  context  the  con¬ 
tribution  of  the  analytical  data  presented  in  earlier  sections  of 
this  report  in  which  physical  principles  were  employed  to  estab¬ 
lish  a  quantitative  relationship  between  free-field  blast  parameters 
and  the  aerodynamic  characteristics  of  objects  that  may  be  dis¬ 
placed  by  blast  winds.  Thus,  one  may  determine  or  estimate  many 
of  the  important  physical  factors,  and  the  quantitative  values  associ¬ 
ated  therewith,  that  are  pertinent  to  the  assessment  of  environ¬ 
mental  hazards. 

For  example,  it  is  desirable  to  know  what  the  velocity  of 
debris  may  be  as  a  function  of  yield,  range,  and  distance  of 
travel  for  inanimate  objects  having  various  areas,  masses,  and 
drag  coefficients.  Likewise,  it  is  of  value  to  know  the  order  of 
magnitude  and  duration  of  the  "G"  loads  imposed  on  animate  objects 
by  blast  winds  associated  with  different  overpressures  produced  by 
various  explosive  yields.  Also,  under  similar  circumstances,  it  is 
helpful  to  have  values  for  the  velocity  of  animate  objects  as  a  function 
of  time  and  distance  of  travel.  The  latter  is  often  pertinent  because 
the  work  space  of  one  exposed  individual  may  allow  only  a  few  feet 
of  travel  and  thus  limit  the  impact  velocity;  for  another  individual 
the  environment  may  allow  attainment  of  a  higher  and  perhaps  maxi¬ 
mal  velocity  before  decelerative  events  occur. 

Thus  the  graphic  data  prepared  for  the  present  study  not  only 
contribute  to  the  physical  aspects  of  blast  effects,  but  also  offer 
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information  of  value  to  thos  interested  in  blast  and  shock  biology 
as  will  be  noted  briefly  below. 

4.  3.  4  Biological  Interests 

There  are  at  least  two  reasons  why  quantitative  data  relevant 
to  blast-induced  translation  of  objects  interest  biomedical  person¬ 
nel.  The  first  is  entirely  pragmatic,  but  requires  that  enough  in¬ 
formation  about  biologic  response  be  available  to  formulate  biologic 
criteria  equal  to  the  challenge  of  hazards  assessment.  When  such 
criteria  exist,  it  becomes  analytically  possible  to  set  forth,  as 
functions  of  yield  and  range,  "safe"  areas  and  those  within  which 
performance  may  be  degraded,  casualties  may  occur,  and  various 
levels  of  lethality  can  be  expected. 

The  second  reason  physical  data  relevant  to  blast-induced 
environmental  variations  intrigue  blast  biologists  is  related  to 
the  fact  that  biological-response  data  are  frequently  lacking  or  are 
inadequate  for  hazards  assessment.  Under  such  circumstances 
the  physical  information  can  be  used  to  plan  conceptually  and  to 
direct  more  realistic  research  programs.  A  case  in  point  concerns 
the  very  high  initial  G-loads  predicted  for  objects  the  size  and  shape 
of  man  set  forth  in  Table  3.  2  and  the  G-time  patterns  contained  in 
Figs.  3.  25  to  3.  36  applicable  to  the  168 -lb  man;  viz.  ,  acceleration 
coefficient  (a)  values  of  .  052,  .  021,  and  .  0063  for  individuals 
standing  facing  the  wind,  crouching  facing  the  wind,  and  prone 
aligned  with  the  wind,  respectively.  The  physical  data  strictly 
refer  to  the  displacement  of  the  center  of  gravity  of  rigid  objects 
simulating  an  "average"  man.  They  specify  G-loadings  that  rise 
"instantaneously"  to  very  high  values  and  decay  differently  with 
time  depending  upon  yield,  range,  acceleration  coefficient,  etc. 

They  say  nothing  about  the  G-time  variations  that  actually  occur 
on  the  down-stream  side  of  semi-elastic  living  object  compared 
with  the  up-stream- side  or  about  the  associated  loads  applied  to 
different  internal  body  organs. 

The  physical  data,  however,  do  pose  problems  for  perceptive 
biologists.  For  example,  what  is  the  biology  of  instantaneously 
applied  G-loadings?  Can  high-density  blast  winds  produce  injury 
only  because  they  suddenly  "push"  a  man  too  fast,  and  if  so,  under 
what  circumstances?  Are  the  significant  effects,  if  any,  limited 
to  small  explosive  charges  and  to  "isolated"  portions  of  the  body 
such  as  fingers,  feet,  extremities,  etc.?  What  is  the  comparative 
range-yield  relationships  between  these  kinds  of  G-loads  and 
hazards  due  to  primary  (pressure)  and  secondary  (missiles)  blast 
effects  ? 


-59- 


L 


